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£H , Abstract. We use Green's function method to study the Heisenberg model of 

LaFeAsO with the striped antiferromagnetic collinear spin structure. In addition to 
the intra-layer spin couplings J\ ai Jib, J2 and the inter-layer coupling J c , we further 
consider the contributions of the single-ion anisotropy J s . The analytical expressions 



U 

Oh 

^ ' for the magnetic phase-transition temperature Tjv and the spin spectrum gap A are 

1 obtained. According to the experimental temperature T/v = 138K and the previous 

estimations of the coupling interactions, we make a further discussion about the 
magnitude and the effects of the single- ion anisotropy J s . We find that the magnitudes 
of J s and J c can compete. The dependences of the transition temperature Tjy, the 
zero-temperature average spin and the spin spectrum gap on the single-ion anisotropy 
are investigated. We find they both increase as J s increases. The spin spectrum gap 
O . at low temperature T — > is calculated as a function of J s , the result of which is a 

useful reference for the future experimental researches. 
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1. Introduction 

It was recently discovered that an iron-based material LaFeAsO shows high-temperature 
superconductivity when atoms are partially substituted by F atoms [I]. This discovery 
has triggered great research interest on the FeAs-based pnictides superconductors and 
their undoped compounds. It has been theoretically and experimentally confirmed 
that these pure FeAs-based compounds have a ground state with collinear stripe-like 
antiferromagnetic(AF) spin order formed by Fe atoms[2l El lU El EJ El EJ [101 E] ■ Thus 
to establish an effective spin Hamiltonian for them and to elucidate the corresponding 
antiferromagnetism are helpful in understanding the underlying mechanism to make 
them superconducting upon doping. 

For undoped LaFeAsO and other similar parent compounds, a Heisenberg exchange 
model was suggested to explain their AF structure [21 Efl H2J [TS], and was used to explore 
their magnetic properties [H f [T5| [16] . Figured] shows the unit cell of the orthorhombic AF 
spin structure of the Fe lattice. This orthorhombic structure exists below a structure 
transition temperature Tg, which is 15 ~ 20K higher than the magnetic transition 
temperature T/v[Sl HZ]- Usually the nearest neighbor (NN) coupling J x (including J la 
and J\b), and the next-nearest neighbor (NNN) coupling J 2 in FeAs layers are dominant 
and must be considered. The NN coupling J c between spins on neighboring layers is 
regarded to be much smaller than the in-plane couplings^ Hj. However, J c was found 
to be essential for the existence of a non-zero magnetic transition temperature r^|15j. 
A further consideration can include the single-ion anisotropy J s . It was estimated to 
be even much smaller than J c in the model of SrFe2As2, but a spin spectrum gap was 
found to be produced by it[14j. For LaFeAsO, so far there is no research report about 
the magnitude or the effects of the single-ion anisotropy. 

In this paper, we use Green's function method[18j to study the Heisenberg model 
of FeAs-based pure parent compounds. The Hamiltonian of this model in a detailed 
form is 

1 1 
H = 2^ la X] ^ u ' ~^ 2^ la ^ ' "~' 2jf ^ ^ u ' "~' 2j 

(ij) " (ij) (ij) 

+ J2 £ Su-S 2j + J c J2Su-S 2j - JsYHi^uf + (S z 2i )% (1) 

({ij)) (ij) i 

where the spin coupling J c between layers and the single-ion anisotropy J s are both 
considered. The subscripts 1 and 2 mean the sublattices 1 and 2 respectively, (ij) 
means NN spin pairs, and ((ij)) means NNN spin pairs. The self-consistent equations 
for the average sublattice spin will be derived. An analytical expression for the magnetic 
transition temperature T/v will be obtained. For LaFeAsO, according to the recent 
estimations of the strengths J ia , Ju, J 2 and J C [IS] with the experimental temperature 
T/v = 138K[6l[l7], we shall make a further estimation of the single-ion anisotropy J s . We 
find that the magnitude of J s can compete with J c , and in some situations even bigger 
than J c . The effects of the single- ion anisotropy on the transition temperature T/v, the 
zero-temperature average spin (S z )o and the spin spectrum gap are investigated. We 
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Figure 1. A unit cell of the orthorhombic Fe spin lattice, a, b, c are the three base 
vectors. The lattice consists of two sublattices (distinguished by the color gray and 
black) . 

find they all increase as J s increases. In section |2], we shall give our analytical results 
derived from the Green's function method. In sectional we shall present our numerical 
results. Finally a conlusion is given in section HI 



2. Green's function derivation 



According to the general scheme of Green's function method to solve an 
antiferromagnetic spin model with two sublattices, we construct the following Green's 
functions: 

G lk>u {u) = «S+ ; S u )), G 2kAl (u) = «S+ ; S u )). (2) 

The equation of motion is 

u((A; S u )) = ([A, S u }) + (([A, H}; S u )), (3) 

where A represents the spin operator S^ k or S^- The commutator [A, H] can be derived 
using Hamiltonian (0Q) and the basic commutation relations of spin operators: [Sf, Sj] = 
2Sf6 lJ} [Sf, Sf] = TS?5 ij} where S? = Sf ± iSf. 

In order to close the system of equations, the so-called PRA or Tyablikov 
decoupling [18] is adopted for the terms stemming from the exchange couplings: 

((S*S+;Sr))*(S*)((S;;Sr)), i^j. (4) 

While for the terms stemming from the single-ion anisotropy, we adopt the Anderson- 
Callen(AC) decoupling [19]: 

((St St + S+S*; Sr)) « 2(^)9^ «S+; S-)), (5) 

where 

®l z) = l-^- 2 [S(S + l)-(S?S?)]. (6) 

The AC decoupling has been demonstrated to be most adequate for the single-ion 
anisotropy much small compared to the exchange interactions [201 EI] • 
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In order to write the decoupled equations of motion in the k space, we take the 
following Fourier transformation: 

1 



Gki(oj) 



N 



uj)e 



(7) 



where iV is the number of sites in either sublattice, and the summation over k is 
restricted to the first Brillouin zone of the sublattice. At the same time, the equation 



h J2k e lk '( Ri R ^ is also used. 



N 



Because of the translation invariant, we have (S' k ) = (S z ), (S'ffc) = ~($ z ) an d 



U 2fc 



L [S(S 



2S 2 



(S Z S Z )}. Finally, we obtain the decoupled 

(8) 

[u + (S z )A k ]G 21 (k,u) + (S z )B k G n (k,u) = 0, 
A k = 2J la cos(k x a) - 2J la + 2 J lb + 4J 2 + 2J C + 2J s Q {z \ 



equations of the two Green's function in k space: 

[u- {S z )A k )G n {k,co) - (S z )B k G 21 (k,u) = 2(S Z ), 

and 



where 



(9) 
(10) 



and 



B k = 2Ju,cos(k y b) + 4J 2 cos(k x a) cos(k y b) + 2 J c cos(k z c), 

in which a, b, c are the three lattice constants. Solving equations (j§J) and 
the Green's function: 

Gn(k,u) 

and the spin spectrum: 

U k =(S*)yfA 



A k (S z )+u k A k (S z )-co k 



u -u k 



CO + u k 



Bl 



(11) 
we obtain 

(12) 
(13) 

A = 2(S z )p s Q( z )[2J lb + 4J 2 + 2J C + J a 8W], (14) 

which is similar with the expression given in ref [Hj derived from the spin-wave theory, 
expect for the factor 0". From this expression for the gap, one see that the single-ion 
anisotropy is essential for the existence of the spectrum gap. 

Then following the process of solving the average spin, we derive the correlation 
function (S~S + ) using the spectrum theorem: 

1 



When k — > 0, we obtain a expression for the spectrum gap: 



(s-s + ) 
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00 , ImGn(fc, 00 + ie) 
duj — 
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in which the equation 
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A\ ~ Bl 



coth ^— - — 1 



(15) 



in8(x)(P(- ■ •) means taking the principle value) 



has been used to obtain the imaginary part of Gn(w + ie), and (3 
Boltzmann constant, T is the temperature. 
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ks is the 
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According to the theory of Callen|22j, the average spin for arbitrary 5 can be 
calculated using the following equation: 

(5 -$)(! + $) 25+1 + (5 + 1 + $)$ 25+1 



(S z ) 



'1 + <&) 2S + 1 - $ 2 5+l 



(16) 



where 



(5~5+) 
2(5 2 > 

-E 



.4, 



coth 



(17) 



On the other hand, the correlation function(5 z 5 2 ) can be calculated from the equation 
(5*5*) = 5(5 + 1) - (1 + 2$)(5 2 ). Using equation p, we can relate 6^ to $ by 



e« = i 



25 2 



2$). 



(18) 



Now the equations (fl6l (ITTj) (1TB1 can be solved self-consistently to obtain the average 
spin at any given temperature, provided we know the values of the exchange couplings 
Jia, Jib, J 2, J c and the single-ion anisotropy J s . 



When the temperature T approaches zero, we obtain coth(^p-) — > 1. The equation 



( 1T7|) is reduced to 



T-s-0 



2N 



E 



,4, 



Bl 



- 1 



(19) 



The zero-temperature average spin (5 z )o can be obtained by self-consistently solving 
the equations (ITBD (ITSl) (TT9l) . 

When the temperature T approaches the magnetic transition temperature Tjv, the 
average spin (5 Z ) as well as the spectrum ujk will approach zero. Expanding coth(^ fe ) 
in the equation (ITT)) , we obtain 

r i 

\ t ^ Tn « - -, (20) 



$1 



/?(5 2 ) 2' 

where T = ^ J2k a^-b 2 • inser ting (12171) into (1TB1 . and expanding the terms in the 
denominator and the numerator as the series of (5 Z ), we finally derive 



(S z ) 



where 



Ti 



N 



12(Tk B T N ) 2 
\\ 5(25-1) 

5(5 + 1) 



T 

TV 



3k B T 

On the other hand, inserting (12171) into (ITS]) , and using the equation 
reduced expression for O^ 2 -* near the temperature Tn- 

25-1 



W 



35 



(21) 
(22) 

we obtain the 
(23) 
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Figure 2. The transition temperature T/v as a function of the single- ion anisotropy 
J s for J la = 0.98, J 2 = 0.52 and J c = 0.0004 in the unit of J lb = 50 mcV. The black 
point correspond to the experimental temperature T/v = 0.238 in the unit of Ju,/kg. 



3. Numerical results and discussions 

So far there is no consensus on the magnitudes of the exchange couplings J\ a , Ju, J 2 and 
J c , because of the unclear microscopic origin of the observed AF spin structure. Here 
we prefer the estimations in ref[15], which gave J\b = 50 ± 10 meV, J\ a = 49 ± 10 meV, 
J2 = 26 ± 5 meV and J c = 0.020 ± 0.015 meV by using the experimental transition 
temperature Tn = 138K of pure LaFeAsO. The main purpose of this paper is to 
investigate the magnitude and the effects of the single-ion anisotropy J s in LaFeAsO. 
Through out our numerical calculation, we take Jib = 50 meV, J\ a = 49 meV, J2 = 26 
meV and the spin S — 1. The result Jib ~ 50meV are obtained from the first-principle 
calculating [U [23]. In the present systems of units, the Boltzmann constant is taken as 
k B = 0.086 meV/K. 

Figure [2] shows the effect of the single-ion anisotropy J s on the transition 
temperature T/v- We can see that Tjv increases as J s increases. This means that the 
single-ion anisotropy term is in favor of the AF spin structure. It can be understood 
from the expression of the single- ion anisotropy term in the Hamiltonian ([T]). Increasing 
the magnitude of J s will make the spins incline to align along the z axis, and give a lower 
total energy, which make the system more stable. To one's surprise, the magnitude of J s 
corresponding to the experimental transition temperature Tjv is about 0.00143 Ju, which 
is much bigger than the magnitude of the exchange coupling J c = 0.0004J lb estimated 
in ref[15]. Furthermore, we see from figure [3] that the variation range of (S z )o with J s 
varying from to 0.1 is almost the same as the one produced by J c in ref[15]. All these 
results imply that the magnitude of J s is probably not much small compared with J c . 
So the estimation of J c maybe need to be adjusted if the single-ion anisotropy term is 
considered. As to the estimations of the other exchange couplings J\ a , Jib and J 2 , we 
think there are still reasonable. 
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Figure 3. The zero-temperature average spin as a function of the single-ion anisotropy 
J s for Ji a = 0.98, J 2 = 0.52 and J c = 0.0004 in the unit of J lb = 50 meV. 
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Figure 4. The competing relation of J c and J s . The curve is depicted by using the 
equation (JUD for T N = 0.238(J lb /k B ), J la = 0.98J lb and J 2 = 0.52J lb . 

Figure H] shows the competing relation of J c and J s when the transition temperature 
is fixed at the experimental value. The increase of J c is accompanied by the decrease 
of J s , and vice vera. From figure HI we can see that the ranges of their corresponding 
variations are at the same magnitude, which implies they probably have the same status 
in the viewpoint of theoretical study. As to revealing the actual magnitudes of the two 
parameters J c and J s , we think it is not enough to use only the experimental transition 
temperature T^. 

Figure |5] shows the effect of the single-ion anisotropy J s on the spectrum gap at 
low temperature. The gap vanishes as J s vanishes, and increases as J s increases. We 
find the effect of J c on the gap is very trivial. The curves for different values of J c 
between (0.0001,0.01) are almost the same, while the single-ion anisotropy affects the 
gap apparently. Considering the gap can be obtained from inelastic neutron-scattering 
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Figure 5. The spin spectrum gaps at low temperature T — > as functions of the 
single-ion anisotropy J s for J la — 0.98Jib, J2 = 0.52Jk, and J c = 0.0001 Jif,(the below 
curve), J c — 0.01Jif,(the above curve). 

experiment |14j. we suggest that the magnitude of the single- ion anisotropy be estimated 
from the future experimental results of the spectrum gap. For example, if J c = 0.0004 Ji b , 
we obtain J s = 0.00143Ji b from the experimental transition temperature Tjy = 138K. 
Then calculating the spectrum gap with J s = 0.00143Ji b , we obtain the magnitude of 
the gap A « 3.4 meV, which can be compared with the future experimental result. 

4. Conclusion 

We use Green's function method to study the Heisenberg model ((TJ) of LaFeAsO with 
the striped AF spin structure as shown in figure [TJ The main purpose of this paper is to 
investigate the magnitude and the effects of the single-ion anisotropy J s . We derive the 
self-consistent equations for the average spin, and obtained the analytical expressions 
for the spin spectrum gap A, and the magnetic transition temperature T N . We find 
that the transition temperature Tjv, the zero-temperature average spin (S' z )o and the 
spin spectrum gap A are all increasing functions of the single- ion anisotropy J s . From 
our numerical results by using T N = 138K and the previous estimations of J la , Ji b , J 2 
and J c in ref[15j, we find that the magnitude of J s is probably not much small compared 
with J c . Because the single-ion anisotropy is essential for the existence of the spin 
spectrum gap, we suggest using the experimental result of the spin spectrum gap to fix 
the magnitude of the single-ion anisotropy J s in the future. 
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